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Perfect porcupines are specially-configured networks of gravitational wave detectors, in the limit 
that the individual detectors and the distances between them are short relative to the gravitational 
wavelengths of interest. They have beautiful properties which make them ideal gravitational wave 
telescopes. I present the most important cases explicitly. For a network of one-arm detectors (like 
"AGIS" 01), the minimal perfect porcupine has 6 detectors, oriented along the 6 diameters of a reg- 
ular icosahedron. For a network of two-arm detectors (like the equal-arm Michelson interferometers 
LIGO/ VIRGO) the minimal perfect porcupine is a certain 5 detector configuration. 
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INTRODUCTION 

In this decade, wc hope and expect that gravitational 
waves will be directly detected for the first time, and 
a new field will be born: gravitational wave astronomy. 
Attention will shift, from individual gravitational wave 
detectors, to neworks of multiple detectors which func- 
tion together as a gravitational wave telescope. (For an 
introduction to gravitational wave detector networks, see 
and references therein.) Here I consider the regime in 
which (by choice or necessity) the individual detectors 
and the distances between them are short relative to the 
gravitational wavelengths of interest, so that the arms of 
the various detectors in the network may be thought of 
as emanating from nearly the same point in space, like 
the fanned quills of a frightened porcupine. In this limit, 
the relative positions of the detectors are unimportant, 
and it is their relative orientations that matter. For spe- 
cial configurations (particular choices for the number of 
detectors and their orientations), such networks exhibit 
beautiful properties that make them, in certain respects, 
ideal astronomical instruments. These "perfect porcu- 
pines" are the subject of this Letter. 

Motivation comes from recent developments in gravita- 
tional wave detection. The ground based LIGO/ VIRGO 
detectors ^ have successfully reached their ( "initial" or 
"stage I" ) design sensitivity and will be upgraded over the 
next several years to a higher ( "advanced" or "stage 11" ) 
sensitivity which will likely allow them to detect gravita- 
tional waves. But the initial LIGO/VIRGO experiments 
are only sensitive to high frequency (30 Hz < / < 10"^ Hz) 
radiation; they are analogous to gamma ray telescopes, 
in that their expected signal consists of occasional brief 
bursts which represent only the high frequency tip of the 
astrophysical gravitational wave spectrum. It would be 
tremendously valuable to have complementary detectors 
(analogous to optical or radio telescopes) to study the 
rich array of sources and physical effects at lower fre- 
quencies (/ < 10 Hz). To this end, several space missions 
have been proposed; but, unfortunately, there is uncer- 
tainty about when and whether the first generation mis- 
sion (LISA optimized near / ^ 10"'^ Hz) will launch. 



and even greater uncertainty surrounding the proposed 
second generation missions (BBO [H and DECIGO 0], 
optimized near / ~ 0.3 Hz). Thus, it makes sense to give 
some consideration to the possibility of pursuing low fre- 
quency (/ ^ 10 Hz) gravitational wave astronomy from 
the ground - as best we can, and hopefully in parallel 
with the space missions mentioned above. But then we 
quickly find ourselves in the "porcupine" regime, because 
frequencies / ^ 50 Hz correspond to wavelengths longer 
than the radius of the Earth! 

Recent developments have invigorated the possibil- 
ity of gravitational wave astronomy at low frequencies 
(/ ^ 10 Hz) from the ground. In particular, there are 
ideas for how to design detectors that are less suscepti- 
ble to the two noise sources (seismic noise and gravity 
gradient noise) that limit LIGO/VIRGO at low frequen- 
cies. One idea is to operate underground, where seismic 
noise and gravity gradient noise are intrinsically weaker 
and. moreover, may be more effectively monitored (with 
seismometers) and removed 01 ■ A second idea is gravita- 
tional wave detection using atom interferometry [H H, Q , 
which may be more immune to seismic noise, and has 
many possible avenues for future development. In the 
Discussion, I suggest another idea which may be relevant 
to gravity gradient noise in certain contexts. 



PROPERTIES OF PERFECT PORCUPINES 

In this section I define perfect porcupines and highlight 
some of their good features. In the next two sections, I 
present the most important perfect porcupines explicitly. 

First let me fix notation. My fourier conventions are: 

/ + 00 /• + 00 

#e-2"/*^(/), 5(/)= / dte+2"/*g(0. (1) 
-oo J —oo 

I work in transverse-traceless (TT) gauge [l^l ■ The lower 
case latin indices {i,j} label the 3 spatial directions: 
i,j = 1,2,3. The upper case latin indices {A,B} la- 
bel the 2 gravitational wave polarizations: A,B =^ 1,2. 
The lower case greek indices {a, /3} label the TV detectors 
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in the network: a, f3 = 1, . . . , N. I use the Einstein sum- 
mation convention: repeated indices (one upper and one 
lower) are to be summed. Hats denote unit vectors. 

When gravitational waves reach us from a distant as- 
tronomical source, they appear as a gravitational plane 
wave travelling in the h direction 



A=1.2 



(2) 



where the two polarization waveforms hA^r) are arbi- 
trary functions of t = t~n-x, and the polarization tensors 
Pf^{n) form an orthonormal basis on the 2-dimensional 
space of symmetric, transverse, tracclcss 3x3 matrices: 



[P'i{n)YP^^{h) 









(3a) 
(3b) 
(3c) 
(3d) 



The output s^it) of detector a has two parts, gravita- 
tional wave signal h^{t) and noise n^it): 



(4) 



Let us first consider h^{t) and the information it con- 
tains. In the porcupine limit, all of the detectors in 
the network live at essentially the same spatial location 
(x ~ 0), and hence only measure hij{t,x) at that point: 
hij {t, 0) . The gravitational wave signal is 



hair) 



+ 00 



dTW^{T)h,j{T^T)A:^ 



(5) 



where the window function Wa(T) describes the tempo- 
ral response of detector a, and the "antenna pattern" A^J 
[normalized as {Aaij)*A^^ = 1] is a projector describing 
its spatial orientation. Since each detector a is sensitive 
to a particular linear combination hij of the compo- 
nents of hij, the collection of detectors in the network 
measure linearly independent combinations if and only if 
the determinant of the N x N grammian matrix of the 
antenna projectors is non-zero: 



DetK^ Af,] ^ 0. 



(6) 



We can have up to 6 linearly independent antenna pro- 
jectors, and in this maximal case the network measures 
all 6 components of hij{t,0). Then, even though these 
measurements are confined to x = 0, they are enough 
to determine the properties of the plane wave as fol- 
lows. First, the propagation direction n is the zero eigen- 
vector of hij(t,0). (Actually, this only determines h 
up to a sign, since — n is also a zero eigenvector [Tlj.) 
Then, by projecting hij{t,0) onto the polarization ten- 
sors P^(h), we obtain the two polarization waveforms: 



hA{t) = [P'2{h)]*h^j{t,6). We started with 6 functions 
of time hij{t, 0); so far we have used these to completely 
determine the properties of the gravitational plane wave 
[two functions hA{t) and two angles to specify n\] and 
the remaining information may now be used to perform 
the following cross checks that we are really observing a 
standard gravitational plane wave: 



dh/dt 
Dct[/i,y(t,0)] 



Ti-[h,j{t,0)] = 0. 



(7a) 
(7b) 
(7c) 



Now let us turn to the noise n^(t), which we model 
as stationary and gaussian, with zero mean, so that it 
is characterized by its correlation function Cq^(T) or, 
equivalently, its spectral density Sapif) = Capif)'- 



(8b) 



Sapif) induces a natural inner product on the space of 
signals (or noise) in the network: 



+ 00 



df9i'^*{f)[S-HfT^9f\f). (9) 



A particular noise fluctuation has probability propor- 
tional to exp[— (n|n)/2] and the expected signal-to- noise 
ratio (SNR) of a gravitational wave is (/i|/i)^/^. If a grav- 
itational wave signal (which depends on various param- 
eters dk) is detected, and the likelihood function may be 
approximated as a gaussian oc exp[—{l/2)6kT''''9i] near 
its peak, then the expected inverse covariance matrix is 



dh 



dh 



(10) 



Consider a network of identical uncorrelated detectors: 
WaiT) = W{T) and Sapif) = S{f)6ai3. Such a network 
is a "perfect porcupine" if {h\h) simplifies to the form 



{h\h) = C 



-00 ^ S{f) 



.A=l 



(11) 



where C is a constant. This says that the network's grav- 
itational wave sensitivity is independent of the direction 
or polarization of the wave. 

Perfect porcupines also have other nice properties. 
First, as we shall see, they: (i) determine the propagation 
direction n and both polarization waveforms hA{t) of a 
gravitational plane wave, as described above; and (ii) per- 
mit the systematic checks (ffal I7bp . In addition, perfect 
porcupines built from one-arm detectors (like AGIS [H) 
also permit the other systematic check l\7c\ : we return 
to this point in the Discussion. Next recall that ha{t) 
depends on two sorts of parameters: (i) the two angles 
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in n, which we denote by the labels fi and v; and (ii) 
all other parameters (such as the masses and spins and 
inclinations in an inspiraling binary black hole) which we 
denote by the label a. With this notation, if we choose 
angular coordinates that are "nice" near the point n, in 
the sense that they run along two perpendicular great 
circles through h (like ordinary polar coordinates and 
(j) at the equator), then we find: 

r^"^ = {h\h)S^"' = SNR^S'"' (12a) 
r^'" = 0. (12b) 

This says that the expected uncertainties in the two an- 
gular coordinates of the source are equal to each other, 
mutually uncorrelated, and independent of n; and also 
that they are uncorrelated with the uncertainties in all 
of the other parameters characterizing the source. Fur- 
thermore, the perfect porcupine's angular resolution SO = 
1/SNR should be compared with the angular resolution 
from triangulation: S9 ~ (A/L)(1/SNR), where A is the 
gravitational wavelength and L is the distance between 
the detectors in the network. Thus, in the (porcupine) 
regime where L is short relative to A, the angular resolu- 
tion of a perfect porcupine is parametrically better than 
the angular resolution from triangulation. 



NETWORKS OF ONE-ARM DETECTORS 



independent in the sense of Eq. ([6]): they measure all 6 
independent components of hij, but do so redundantly. 
All else being equal, this redundancy makes the 10-arm 
perfect porcupine better (but also more expensive) than 
its 6-arm counterpart: the redundancy allows us to cross 
correlate two noisy data streams which contain the same 
gravitational wave signal, but different and uncorrelated 
noise, and thereby extract the true signal better than we 
could from either data stream individually. 

An even larger and more redundant perfect porcupine 
consists of 15 detectors, oriented along the 15 directions 
connecting opposite edges of a regular icosahedron (or, 
equivalently, the 15 directions connecting opposite edges 
of a regular dodecahedron). This has C = 2. 

If one relaxes the requirement that all of the detectors 
be identical, more options become available. For exam- 
ple, for > 5, the configuration 



{0,0,1} 



cos^,sin^,i} (a 



(a = 0) 

l,...,fc) 



with 



.JkJZW{T) (a = 0) 
W{T) (a = l,...,fc). 



is a perfect porcupine with C = 4fc/25. The k 
is just the minimal perfect porcupine ([T4| . 



(16) 



(17) 



5 case 



One-arm detectors (such as AGIS [l|) have antenna 
projectors of the form: 



4" 



(13) 



The minimal perfect porcupine built from such detec- 
tors has 6 arms, oriented along the 6 directions connect- 
ing opposite vertices of a regular icosahedron (or, equiv- 
alently, the 6 directions connecting opposite faces of a 
regular dodecahedron): 



{0,0,1} 
\/475{c 



, sm- 



i} {a. 



:0) 

a,. 



.,5) 



(14) 



This network has C ~ 4/5. These 6 detectors are inde- 
pendent in the sense of ([51) ; so in addition to determining 
the direction of a gravitational plane wave n, and both 
polarization waveforms Ha (t) , they also permit the three 
cross checks (|7al I7bl [7c|) to be performed. 

A larger perfect porcupine has 10 detectors, oriented 
along the 10 directions connecting opposite vertices of a 
regular dodecahedron (or, equivalently, the 10 directions 
connecting opposite faces of an icosahedron): 



2ip 



±1/2 



27ra . 27ra 1 ^ ip^'^ 

cos , sm , 

5 ' 5 ' 2 



(15) 



where a = 1, . . . , 5, ip = {1 + ■\/5)/2 is the golden ratio, 
and C = 4/3. Of course, these 10 detectors are not 



NETWORKS OF TWO-ARM DETECTORS 

Two-arm detectors (Michelson interferometers like 
LIGO/VIRGO, with equal and orthogonal arms) have 
antenna projectors of the form 



ifit) 



If we define k = iarccos-\/3/5 and 



a" 


= { 




= { 


then, for a - 


- I,-. 




pa 







^ 27ra 



-cosK a" -f- siuK 6" 
-sinK a" -I- cosK 6" 



(18) 



(19) 



(20) 



is a perfect porcupine with C ~ N/5. This network mea- 
sures the 5 traceless components of hij (with greater re- 
dundancy when N is larger) . It therefore determines the 
direction n of a gravitational plane wave, and both po- 
larization waveforms hAit); and it permits the first two 
cross checks ([2al[7b|, but not the third When N = 5 
this is the minimal perfect porcupine. 

Another perfect porcupine is worth mentioning. Con- 
sider the 15 directions connecting opposite edges of a 
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regular icosahedron (or, equivalent ly, the 15 directions 
connecting opposite edges of a regular dodecahedron). 
These 15 directions separate into 5 orthonormal triads. 
From each orthonormal triad, we can select 3 different 
orthonormal pairs {pa, Qa}- In this way, we obtain a per- 
fect porcupine with 15 detectors, and C = 3. 

DISCUSSION 

We have seen that a perfect porcupine built from one- 
arm detectors will monitor Tr[/iij]. As mentioned above, 
the vanishing of this channel is a check that one is observ- 
ing standard gravitational waves, as opposed to noise, or 
something more exotic. Alternatives to general relativity 
often give the graviton a zero-helicity component; and 
zero-helicity gravitational waves would show up as fluc- 
tuations in Tr[/iij]. Also, as one goes underground, and 
to lower frequencies, the coherence length of the fluctu- 
ations in the Newtonian gravitational potential becomes 
longer: if the various detectors in a perfect porcupine 
can be placed sufficiently close together, so that they all 
see the same (or similar) fluctuations in the Newtonian 
potential, then these fluctuations will also look like fluc- 
tuations in Tr[/iij], and monitoring this channel may even 
be helpful in subtracting gravity gradient noise. On the 
other hand, if it is necessary, e.g., to build the detectors 
parallel to Earth's local gravitational field (so that the 
various detectors in the porcupine must be located at 
widely separated points on the Earth) , then this method 
for removing gravity gradient noise won't work. 

I have focused on a perfect porcupine's ability to mea- 
sure a single plane wave; but this analysis also applies to a 
sum of many plane waves, as long as they are separable in 
the time-frequency (or template) domain. For example, 
two plane waves of different frequencies are not a problem 
(they may be cleanly separated in frequency space, and 
then handled independently); but two plane waves mov- 
ing in different directions at the same frequency would 
be a problem: the cross check (j7bp and the porcupine's 
direction-finding algorithm would fail. Fortunately, al- 
though source "blending" (non-separability) docs occur 
for LISA sources when / < 10~^ Hz, it is unlikely that a 
ground-based perfect porcupine would reach sufficiently 
low frequencies and sufficiently good sensitivities for this 
to be a practical concern. 



An earth-bound porcupine must contend with (in or- 
der of importance) the spin of the Earth, the Moon's 
orbit around the Earth, and the Earth's orbit around 
the Sun; but above we considered perfect porcupines 
which move in straight line, without rotation. As long 
as we are looking at gravitational waves with periods 
much shorter than a day (which is the realistic case), 
then this is a good starting point, for the same reason 
that for many purposes LIGO/VIRGO may be modeled 
as moving through flat space: over timescales contain- 
ing many wave cycles, the detectors are moving along a 
nearly straight non-rotating trajectory. Nevertheless, it 
is important to extend the above analysis to include the 
rotation of the Earth, etc., especially for the purposes 
of studying sources that are detectable over timescales 
longer than a day. It will also probably be best to align 
one of the perfect porcupine's symmetry axes with the 
Earth's rotation axis, if possible. 

Finally, the considerations in this paper were moti- 
vated by recent developments (outlined in the Introduc- 
tion), but I should also mention a more futuristic possi- 
bility: we may eventually be led to build low- frequency 
gravitational wave detectors on the Moon, where the seis- 
mic and gravity gradient noise levels are much lower than 
on Earth. This would again lead us to porcupines. 

I acknowledge support from the CIFAR JFA. 
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